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Abstract. We discuss properties of root posets for finite crystallographic root 
systems, and show that these properties uniquely determine root posets for the 
noncrystallographic dihedral types and type H3 , while proving that there does 
not exist a poset satisfying all of the properties in type H4. We do this by 
exhaustive computer searches for posets having these properties. We further 
give a realization of the poset of type H3 as restricted roots of type Dg, and 
conjecture a Hilbert polynomial for the q, t-Catalan numbers for type H4. 
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1. Introduction 

Let $ be a finite root system with simple roots A and positive roots <I> + . For 
crystallographic root systems, the root poset is given by the partial order on positive 
roots defined as 

/3</3' if /3'-/3eNA. 
Multiple enumerative properties of this poset and in particular about the collection 
of antichains in this poset can be deduced from the degrees of $, which are defined 
to be the degrees of the fundamental invariants for the Weyl group W($) associ- 
ated to <£>. We collect properties of root posets below in Properties 1 through 6. 



Most of the results of this article were achieved at the Leibniz Universitat Hannover in summer 
2012. 
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Even though all the enumerative information can be studied as well for noncrystal- 
lographic root systems, the analogous posets in these cases with N in the definition 
being replaced by K>o do not have the desired properties. 

This led D. Armstrong to raise the question whether there exist posets for the 
noncrystallographic types having these properties. He moreover positively answered 
the question for the noncrystallographic dihedral types and for type H3, see [Arm06, 
Section 5.4.1]. Even though he did not take Properties 5 and 6 into account, the 
posets he constructed do indeed fulfill also those. 

This leaves type H4 as the only type remaining. The main purpose of this paper 
is to discuss this type in full generality. 

Theorem 1.1. Possible root posets for noncrystallographic types are given as fol- 
lows: 

(i) The posets in Figure 1 on page 8 are the unique posets in types ^(m) and H% 

satisfying Properties 1 through 6. 
(ii) There are billions of posets in type H4 satisfying Properties 1 through 4- 
(Hi) None of these posets in type satisfy Property 5. 
(iv) None of these posets in type H± satisfy Property 6. 

Remark 1.2. We want to remark that even in the crystallographic types, Property 4, 
Property 5(b), and Property 6 are partially conjectured and remain open in full 
generality. 

Remark 1.3. In [CK07], Y. Chen and C. Kriloff studied properties of root posets 
in noncrystallographic types in the sense of replacing N by K>o in the definition as 
mentioned above. Their study has a somewhat different flavor since they study the 
correspondence between antichains in root posets and region in Shi arrangements. 
As also mentioned above, these posets do not have the desired properties in the 
sense we discuss here. Nevertheless, the authors obtain, similarly to our results, 
well behaved connections between antichain and these regions in types Izipi) and 
H3, while they show that the connections in type if 4 are not well behaved. 

We prove Theorem 1.1 by an exhaustive computer search for posets satisfying 
the various properties. We used (and present here) two different algorithms for 
the computations. In both algorithms we construct posets inductively. The break- 
through is achieved by testing the right properties at the right moment. 

In Section 2, we review the combinatorial properties of root posets for finite 
crystallographic root systems. In Section 3, we prove Theorem l.l(i). Section 4 
is devoted to provide two algorithmic approaches to prove Theorem 1.1 (ii)— (iv). 
In the final Section 5, we describe how we obtained the conjectured Hilbert series 
needed to provide Property 6 for type H4. 

2. Root posets and their combinatorics 

Throughout this paper, let $ be an irreducible and finite root system of rank n 
with simple roots A and positive roots $ + . Moreover, let W = W($) be the 
corresponding reflection group acting on a real vector space V, c £ W be a Coxeter 
element, d\ > ■ ■ ■ > d n be the degrees for W, and h = d% its Coxeter number. We 
refer to [Hum90] for background on finite root systems and reflection groups. It is 
well known that irreducible finite root systems are classified according to Cartan- 
Killing types 
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• A n with n > 1, 

• B n with n > 2, 

• C„ with n > 3, 

• D n with n > 4, 

• G2,F4,E 6 ,Ei,E s , 

• #3, -H4, and ^(iti) with to = 5 or m > 7. 

A root system $ is crystallographic if the group W($) stabilizes a lattice in V. It is 
well known that the only noncrystallographic hnite root systems are those of types 
H and I. 

For a crystallographic root system, define the root poset $ + = ($ + ,<) to be 
given by the partial order j3 < /3' :<^> /?' — (3 £ NA on positive roots. This poset is 
graded, i.e. all maximal chains in $ + have the same length. The rank of a positive 
root /3 = J^aGA ^aOi is given by X)q<ea ana - ^ * s known that the maximal rank 
is equal to ft — 1. The following definitions make sense for any poset, though we 
discuss them here only in the context of root posets. For a set X C $+, we denote 
by min(X) and by max(A) the set of all minimal and maximal positive roots in X, 
respectively. An antichain in $ + is a set of positive roots of mutually incomparable 
elements. Wc denote the collection of all antichains in $ + by 

A($) = {AC$ + : J 4 antichain}. 

Moreover, an order ideal in <I> + is a set J of positive roots such that 

e I, 7 < P =>■ 7^, 

and we define as well to be the collection of all order ideals in <3? + ($). Obvi- 
ously, antichains and order ideals are in one-to-one-correspondence by sending 

A 1 y 1(A) := {7 G $+ : 7 < /3 for some /3 e A}, 
/(-)• A(J) :=max(J). 

We are now in the position to discuss the enumerative properties of antichains 
(or order ideals) in root posets for crystallographic root systems. In Property 1, 
we collect some direct consequences of the definition, Property 2 was discussed 
e.g. in [Hum90, Theorem 3.20], Properties 3 and 4 can be found in [Arm06, Sec- 
tion 5.4.1], while Properties 5 and 6 are taken from [AST10] and from [StuIO], 
respectively. All these enumerative properties can as well be studied for noncrys- 
tallographic root systems, and each of them can be used to find possible root posets 
in these types. 

Property 1 (Basic properties of root systems). This first property is rather a 
collection of "natural" properties of root posets that are satisfied by definition. 

1. The minimal elements in $ + correspond to the simple system A for the root 
system $. 

2. There is a unique maximal element in the poset ( I >+ . 

3. The poset $ + is graded. 

4. The poset $ + behaves well with respect to standard parabolic subsystems. This 
is, the subposet of $ + given by all j3 € $ + which lie in a subspace generated by 
a subset A' of the simple system A is equal to the root poset for the root system 
<!>' generated by A'. 
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rk($+) 


\A{*+)\ 


him) 


(2, l m ~ 2 ) 


TO + 2 


H 3 


(3,2 4 ,1 4 ) 


32 


H 4 


(4,3 10 ,2 8 ,1 10 ) 


280 



Table 1. The desired ranks and desired numbers of antichains. 



Property 2 (The degrees of a root system). The second property describes how 
the degrees of a root system are encoded in the root poset, or rather in the rank 
generating function of the positive roots. It is taken from [Hum90, Theorem 3.20]. 

Theorem 2.1. Let d\ > . . . > d n denote the degrees o/$. Then the number rkj of 
positive roots in <£> + of rank i is given by 

rk; = {j : dj > i} . 

In particular, the maximal rank in <I> + is given by h — 1 . 

The ranks rk(<I> + ) = (rki, . . . , rk/j_i) of a noncrystallographic root poset $ + 
satisfying this property are collected in Table 1 , where we write i 3 for j consecutive 
ranks of size i. 

Property 3 (Counting antichains). The ^-Catalan numbers are defined as 

di + h 



Cat($) = n ■ 



<=i 

The following theorem is due to A. Postnikov, see [Rci97, Remark 2]. 

Theorem 2.2. Antichains (or equivalently order ideals) in the root poset are 
counted by the $- Catalan numbers, 

\A(&)\ = Cat($). 

This theorem yields the desired numbers of antichain for root posets of noncrys- 
tallographic types as collected in Table 1. 

The following two properties (Properties 4 and 5) are both strictly stronger than 
this property. This is to say that both provide in particular the above counting 
formula for antichains. 

Property 4 (Counting antichains according to size). This property concerns the 
number of antichains of a particular size containing a particular number of simple 
roots. These numbers are encoded in the following generating function 



E 



There is a remarkable conjectured connection between this and two other gener- 
ating functions. These encode information about the noncrossing partition lattice 
and about the cluster complex. The three generating functions were defined by 
F. Chapoton in [Cha04, Cha06] and are called Chapoton's H-triangle, M-triangle, 
and F-triangle. We restrict our attention here to the H- and the Af-triangle. For a 
detailed treatment of all three, see Chapoton's original papers [Cha04, Cha06] and 
as well [Arm06, Section 5]. 
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H($;s,t) 




l + 2st + (m- 2)t + 


s 2 t 2 


H 3 


1 + 3st + 12t + 3s 2 t 2 


+ Ast 2 + 8t 2 + s 3 t 3 


Hi 


1 + 4st + 56t + 6s 2 t 2 


+ 19si 2 + 133t 2 + 4s 3 t 3 + 5s 2 t 3 + 9st 3 + 42i 3 + sH A 



Table 2. The desired 7J-triangles. 



The M-triangle encodes information about the Mobius function on the noncross- 
ing partition lattice NC(&) given by all elements in the interval [1, c] in W endowed 
with the absolute order where c denotes a Coxeter element in W . It is defined by 

Af ($; x, y) = Y, Mfa 7r)x n ~ lk ^ y»- rk « , 

cr,7reA r C(*) 

tT<7T 

where fi : ATC($) x M7($) — > Z denotes the Mobius function on the lattice M7($). 
Observe that we dropped the Coxeter element c from the notations. This is due 
to the fact that the resulting lattice is independent of the chosen Coxeter element. 
F. Chapoton conjectured deep numerical interactions between antichains in the root 
poset and the Mobius function on the noncrossing partition lattice. This conjecture 
remains open. 

Conjecture 2.3. Chapoton's H -triangle is related to the M-triangle via the identity 

s,t) = (1 + (s - l)t T M (*; 7 i T> • 

Conjecture 2.3 can now be used to determine the sizes of the antichains in non- 
crystallographic root posets satisfying this property. The corresponding generating 
functions are collected in Table 2. They are taken from [Arm06, Figure 5.14] and 
independently verified using both the M- and the _F-triangle. 

Property 5 (Counting antichains according to Panyushev orbits). This property 
concerns a cyclic action on antichains in the root poset defined by D.I. Panyushev 
in [Pan09]. He conjectured multiple numerical results which were then later proven 
by D. Armstrong, H. Thomas, and the second author in [AST10]. The Panyushev 
map P : A(Q) — > A($>) is given by mapping an antichain A to the collection of 
minimal elements in the root poset among all elements that are not in the order 
ideal generated by A, 

P(A) :=min (<f> + \ 1(A) ). 

It is shown in [AST10] that the orbit structure of P on A(<& + ) coincides with an- 
other remarkable action on the noncrossing partition lattice given by the Kreweras 
map K : M7($) — ► M7($) which is defined by 

K(cr) = (j- x c. 

This was done by constructing an equi variant bijection ^ : A(§)— —}NC{$). We 
denote by C($) the multiset of sizes of Kreweras orbits in AC($), 

£>($) = {{|0| : O C W is a Kreweras orbit } . 

The multisets of Kreweras orbit sizes for the noncrystallographic types are collected 
in Table 3. They should coincide with the Panyushev orbit sizes of a noncrystallo- 
graphic root poset satisfying this property. 
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h(m) 


I2,m} 




{{2,10 3 }} 


Hi 


{{2,3,5,30 9 }} 



Table 3. The desired multisets of Panyushev orbit sizes. 



The following extensions of this property complete the collection of (partially 
conjectured) properties of root posets. 

Property 5(a). As conjectured in [Pan09] and proven in [AST10], the average 
number of positive roots in an antichain in a Panyushev orbit is constant, 

|o| ^ |A| = ?' 

1 1 Aeo 

where O is any Panyushev orbit in A(<&), and n is the rank of $. 

Property 5(b). The previous property has a counterpart on the average number 
of positive roots in an antichain in a "restricted Panyushev orbit" . Let <&+ s be 
the poset obtained from the root poset $ + by deleting the simple roots, and let 
P res : A($~^, s ) — > A(^^, s ) be the restricted Panyushev map defined in the same 
way as the Panyushev map P : A($ + ) — > A(§ + ) with the poset <3 )+ replaced by 
$+ 5 . It is conjectured in [Pan09] that the average number of positive roots in an 
antichain in a restricted Panyushev orbit is as well constant, 

\0\^ Q l 1 2(^-1)' 

where O is any restricted Panyushev orbit in _A(<I>+ S ), n is again the rank of $, 
and where h is the Coxeter number. Observe here that this conjecture holds true 
for simple reasons that can be deduced from [AST10] in types A and B, and was 
checked in the exceptional crystallographic types. Thus, only type D remains open, 
and we expect that this type can as well be deduced from the results in [AST 10]. 

Property 6 (Counting antichains according to order ideal sizes). The last - again 
conjectured - property comes from the theory of q, t- Catalan numbers. Those were 
extensively studied in type A by various authors. We refer to [GH96, Hag08] and 
the references therein for detailed background and definitions. 

Let W = W{<&) act diagonally on C[x,y] = C[xx, . . . , x n , yx, . . .,y n ] = C[V®V], 
let J be the ideal in C[x, y] generated by all polynomials / g C[x, y] such that 
w{f) = det(ui)/ for all w S W , and define the diagonal coinvariant ring as 

M($) = j/(x,y)J. 

The ($; q 7 t)-Catalan numbers are defined to be the bigraded Hilbert series of the 
diagonal coinvariants (see Section 5 for details), 

Cat($;g,t) =H(M($);q,t). 

In the case of the symmetric group, it was proven by M. Haiman in [Hai02] that 
the specialization t = 1 in Cat(A„; q, t) equals the size generating function of order 
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Cat($;g,i) 


h{m) 


[m + 1] + qt 


H 3 


[16] + qt[W] + qt[Q] 


Hi 


[61] + qt[49] + qt[41] + g 2 f 2 [37] + qt[31]+ 
q 3 t 3 [25] + qH 2 [21] + [13] + q 6 t 6 + q 10 t W 



Table 4. The (conjectured) q, t-Catalan numbers for the noncrys- 
tallographic root systems. 



ideals in the root poset of type A n , 

Cat(A n ;q,l) = ]T gl'L 
iei{A n ) 

In [StulO], the second author studied the ($; q, i)-Catalan numbers for other re- 
flection groups, and conjectured that this phenomenon describes the situation in 
general. 

Conjecture 2.4. The specialization t — 1 in the ($;g, t)- Catalan numbers yields 
the size generating function of order ideals in the root poset, 

Cat($+;<7, 1) = <7 m - 

In Table 4, we provide the (conjectured) ($; q, t)-Catalan numbers for the non- 
crystallographic types, where we write 

[ n ] = [n] q t := ^— — = q^ 1 + q n ~ 2 t +... + qt*-* + t n ~ l . 
q-t 

For the dihedral types him) and type H3, these Hilbcrt series can be directly 
computed using the definition, with the total degree of the polynomials / g C[x, y] 
being bounded by the number of positive roots. The computations for type H4 are 
again much more complicated since the brute force approach is currently far beyond 
the feasibility of modern computers. We discuss these computations in Section 5. 

3. The poset of dihedral types and type H 3 

3.1. The dihedral types. In type I^m), one can easily check by hand that the 
poset shown in Figure 1 (left) is the unique poset satisfying Properties 1 and 2. It 
as well satisfies Properties 3 through 6. 

3.2. Type H3. Let W be the Coxeter group of type H3, i.e. the group with Coxeter 
graph: 

5 

• • • 

We use Algorithm 4.5 below to prove the uniqueness of a graded poset (r = 
{1, . . . , 15}, -<) satisfying Properties 1 through 6. 

Indeed, we have more in this type. Namely, one can construct 1 the root poset 
of type H3 in a natural way from the root poset of type D 6 . Let r = be 



Hugh Thomas informed us that he has a similar approach to the one we use here, in which 
one can construct the root poset of type H3 as a subposet of the root poset of type D§ . 
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Figure 1 . The (unique) root posets of types I2 (m) and H3 



(0,0,0,0,1,1), 

(0,0,0,1,-1,0), 

(0,0,0,1,0,1), 

(0,0,1,-1,0,0), 

(0,0,1,0,0,-1), 

(0,0,1,0,1,0), 

(0,1,0,0,-1,0), 

(0,1,0,0,0,-1), 

(0,1,0,1,0,0), 

(0,1,1,0,0,0), 

(1,0,0,0,0,1), 

(1,0,0,0,1,0), 

(1,0,0,1,0,0), 

(1,0,1,0,0,0), 

(1,1,0,0,0,0), 



(1,-1,0,0,0,0 
(0,1,-1,0,0,0 
(1,0,-1,0,0,0 
(0,0,0,0,1,-1 
(0,1,0,-1,0,0 
(1,0,0,-1,0,0 
(0,0,0,1,0,-1 
(0,0,1,0,-1,0 
(1,0,0,0,-1,0 
(1,0,0,0,0,-1 
(0,0,0,1,1,0) 
(0,0,1,0,0,1) 
(0,1,0,0,0,1) 
(0,1,0,0,1,0) 
(0,0,1,1,0,0) 



Figure 2. The positive roots of the root system of type Dq 



the golden ratio. Let be the root poset for the root system corresponding 
to the Weyl group of type Dq acting on the Euclidean space K 6 with standard 
basis Vi,...,ve and bilinear form (vi \ Vj) — Sij, see [Kac90, §6.7]. Then a\ — 
Vi — V2, • ■ ■ , ot5 — V5 — ve, «6 = t'5 + vq are the simple roots in . With respect 
to the basis v\, . . . ,v&, the positive roots of $J are listed in Figure 2, where a pair 
(a, f3) in a row is given as follows. Write v\ = Y]j aijOLj for ajj € Z. With 

e : Z 6 — > Z 6 , Vi H> (fflj,!, . . . , difi), (base change) 

7 : 1? -)■ Z[r] 3 , (a 1 ,a2,a 3 ,a4,a 5 ,a 6 ) h-» (a x +a 2 r, a 3 + a 4 r, a 5 +a 6 r), 
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Figure 3. Expected number of antichains depending on the num- 
ber of simple roots they contain for root posets of type H4. 

we get T"f{e{a)) = 7(e(/3)). Let <J>^ 3 be the set of positive roots of the root system 
of type H3. With respect to any simple system A, the coordinates of the elements 
of <f> + are in Z[r], Now £~ 1 (7~ 1 ( < 1 > ^.,)) is just the left column of Figure 2. 

In other words, there is a subspace U < K 6 such that the restriction of the 
reflection arrangement of type D e to U is the Coxeter arrangement of type H 3 . 
Under this restriction, pairs of roots in <J>^ g are mapped to pairs a, ra € M 3 where a 
is a positive root of type H3. For each such pair, choose the lexicographically greater 
element in Z 6 with respect to v\ , . . . , v§ . This way, for each element of 4>^ 3 we get 
a choice 1 or r, i.e. a map a : — > {l,r}. The poset ({a(a)a | a € $^ }, <), 
where /?</?' is given by f3' — (3 e N 3 , is (isomorphic to) the poset in Figure 1 (right). 

4. Posets of type ff 4 

Let W be the Coxeter group of type H4, i.e. the group with Coxeter graph: 

5 

• • • • 

Throughout this section, assume that 

(r = {i,...,60},-o 

is a poset satisfying Properties 1 through 4. This is to say that it has 4 minimal 
elements, a unique maximal element, and behaves well with respect to standard 
parabolic subsystems. Moreover, its grading is given as shown in Figure 4, compare 
Table 1. Its number nk, m of antichains A C T of T with m elements such that 
I A n {1, 2, 3, 4} I = fc is shown by Figure 3, where the entry in row k and column m 
is nh,m, compare Table 2. 

4.1. Some consequences of the assumptions. 

Proposition 4.1. We may assume without loss of generality that T is build upon 
Figure 4. I.e., T contains at least the continuously drawn covers and does not 
contain the waved ones. 

Proof. Most of the covers in Figure 4 follow the assumption that T behaves well 
with respect to standard parabolic subsystems. Here, we assume the ordering of 
the simple roots in T to be the same as shown in the above Coxeter graph. 

The following covers have another explanation. The covers in ranks 19 to 29 are 
necessary to ensure the uniqueness of the minimal and maximal elements. For ranks 
12 to 19, we can assume for the same reason that, up to reodering of the vertices in 
the upper rank of two consecutive ranks, that we have at least the following covers 
for each level: 
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Figure 4. Required and forbidden covers in root posets of type 
H4 satisfying Properties 1 through 4. 
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I I 

Again for the same reason, we have the covers 10 -< 13 -< 16 -< 18, and 16 -< 19, 
and as well 18 -< 21 -< 24 -< 27 and 19 -< 22 -< 25 -< 28. 

There has to be a cover for 26 and one for 27. If these are equal, say 29, and if 
there is no further such cover, 26, 27 30, 31, then {26, 27, 30, 31} is an antichain 
with 4 elements, contradicting the assumption that the set of minimal elements 
is the unique antichain of size 4, compare Figure 3. Hence up to permutation of 
29, 30, 31 we may assume 26 -< 29 and 27 -< 30. 

Now we come to ranks 10 and 11: we have a cover for each of 29, 30, 31 and each 
of 32, 33, 34 is a cover. Up to permutation of 32, 33, 34, we have at least 

III- NNL 

But to avoid again an antichain with 4 elements in the second case, we need to add 
a cover. Adding one further cover always yields a situation which may be permuted 
to include the first case, except for 




which again has an antichain with 4 elements. Thus we may assume the covers 
29 -< 32, 30 -< 33, 31 -< 34. □ 

4.2. Assuming as well Property 5. In this section, we moreover assume that 
r satisfies Property 5. First, observe that we can predict two further relations in 
Figure 4, namely the covers 9 -< 13 and the noncover 8 7^ 13. Assume that 8 -< 13. 
Then under the Panyushev action we have the orbit 

{8}, {3, 4, 11}, {6, 7}, {1, 10}, {5}, {3, 4, 8}, {6, 7, 11}, {9, 10}, 

independently of the choice whether we have the cover 9 -< 13 or not. Since it has 
length 8, this would contradict Property 5. On the other hand, removing the vertex 
1 should result in a poset of type A3, which is now only possible if 9 — < 13. 

We present two algorithms to prove the nonexistence of a poset satisfying Prop- 
erties 1 through 4 and also Property 5. Both variants need functions to count or 
to compute antichains of a given poset, and a function to compute the lengths of 
the orbits under the Panyushev action. 

Remark 4.2. We did not use computer algebra systems or similar scripting lan- 
guages for several reasons. Our computations do not require higher functions or 
libraries, and the occurring arithmetic is always with small numbers. Further, since 
our posets have at most 60 elements, we can for instance store the information of an 
antichain or of an order ideal in a very compact 64-bit variable. Last but not least, 
a compiled and optimized version always performs much faster than interpreted 
code. This can be decisive if the runtime is in the range of a few hours to a few 
days. 
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4.2.1. The first algorithm. There are 37 potential covers i -< j with i,j < 36 that 
may be added to Figure 4 together with the cover 9 -< 13 and the noncover 8 -< 13. 
Let us call them ci, . . . , C37. Then there are 14 further potential covers in the upper 
part, say w%, . . . , 1014. We divide the algorithm into two parts. 

The upper part. First we compute all 2 14 possible posets on the vertices 35, . . . , 60 
and store those which do not contain an entire orbit under the Panyushev action 
of length other than 2,3,5, or 30 in a list R. For each element Tq of R, we also 
compute the number O2(ro) of antichains with two elements, the number 335^0) 
of elements not greater or equal than 35, and the number g^e (Po) of elements not 
greater or equal than 36. 

The lower part. For the covers between vertices in 1, . . . , 36, we use the following 
recursion. 

Algorithm 4.3. EnumeratePosetsl(r',p) 

Enumerates all root posets for H4 satisfying Properties 1 through 5, starting from V . 
Input: A poset T' , an index p for the next potential cover. 
Output: Posets based upon V satisfying Properties 1 through 5. 

1. If p = 38 then: 

(1) Count the antichains with 3 elements within 1, . . . , 36. If this number 
is greater than 60, then return 0. 

(2) Compute all orbits under the Panyushev action which are entirely 
contained in 1, . . . , 36. As soon as one of them has length other than 
2,3,5, or 30, return 0. 

(3) Within 5, . . . , 36: Count the number b 2 of antichains with 2 elements, 
the number Z35 of elements not less or equal to 35, and the number 

of elements not less or equal to 36. 

(4) For each T„ e R with a 2 (T ) + b 2 + l 35 ■ g 36 (T ) + l 36 ■ g 35 (T ) = 133, 
compute the poset T as a combination of P' and Po . If P satisfies all 5 
properties, then print T to an output hie. 

2. If all chosen covers between degree 4 and 5 have been included, then check 
that {14, 15, 18, 19} is not an antichain with 4 elements and return if it 
is. 

3. If all chosen covers between degree 8 and 9 have been included, then check 
that 27 -< 29 or 28 -< 29 and that 27 -<! 31 or 28 -< 31 and return if 
not. Further check that {26, . . . , 31} does not contain an antichain with 4 
elements. 

4. If all chosen covers between degree 10 and 11 have been included, then 
check that each of 35, 36 is a cover and that each of 32, 33, 34 is covered. 
Otherwise return 0. 

5. Call EmimeratePosetsl(r',p + 1). 

6. Include the cover c p to V to a new poset T" . 

7. Call EnumeratePosetsl(r",p + 1). 

We use the data of Figure 4 to compute an initial poset V for the first call 
EnumeratePosetsl(r',l). 

Remark 4.4. It is easy to parallelize Algorithm 4.3. For instance, call ci, . . . , cq the 
potential covers within 32, . . . , 36. Now one can fork (duplicate the program) for 
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each of the 25 possible choices for ci , . . . , c% (some of the 2 6 choices are excluded 
by Step 4). 

4.2.2. Results. Algorithm 4.3 terminates after approximately two days of CPU- 
time and does not find any poset satisfying all the above conditions. However, we 
slightly modified the algorithm, and found over 80 billion posets only satisfying 
Properties 1 through 4. We guess that there are several 100 billions of them. We 
show four examples of such posets in Figure 5. Property 5 is much more restrictive 
and turns out to be the best break condition for the algorithm. 

4.2.3. The second algorithm. Since it is very difficult to ensure that an implemen- 
tation of the above algorithm is free of mistakes, we propose in this section a second 
slightly different algorithm, and thereby increasing our changes of providing valid 
results. Further, this second algorithm is much better in the sense that it takes only 
a few hours of CPU-time and performs a more general search, namely for posets 
satisfying all the properties as in the previous sections, except that we do not take 
the structure of standard parabolic subposets into account. For each pair of con- 
secutive ranks we get a set of possible configurations of covers. For example, there 
are 51 possible configurations of covers between two ranks containing 3 vertices, up 
to permutations of the vertices of larger rank: 

iiiNimMiiNmiMmMimmmmmmm^ 

Here, we also assumed that the collection of minimal elements is the unique an- 
tichain containing 4 elements. Similarly, there are 13 configurations for the covers 
between rank 11 and 12, 4 configurations for the covers in ranks 12 to 19, and only 
1 configuration for the ranks 19 to 29. 

Let d be the highest degree, for instance d — 29 in type H4. 

Algorithm 4.5. EnumeratePosets2(r',p) 

Enumerates all root posets starting from T' . 

Input: A poset V , an index p for the next degree. 

Output: Posets based upon T' satisfying Properties 1 through 5, except for the 
standard parabolic substructure. 

1. If p = d then: 

(1) Count the antichains with 2 or 3 elements which do not contain a 
simple root. If this number wrong, return 0. 

(2) Compute orbits under the Panyushev action. As soon as one of them 
has a forbidden length, return 0. 

(3) Print T to an output file, return 1. 

2. If all elements of the last added degree are greater than all simple roots, 
then count the antichains with 2 or 3 elements which contain simple roots. 
These numbers will not change anymore, so we may return if they are 
wrong. 

3. Count antichains and use Property 4 as upper bounds, return if one of 
them is violated. 

4. Compute orbits under the Panyushev action. As soon as one of them has 
a forbidden length, return 0. 
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Figure 5. All root posets for H4 satisfying the Properties 1 
through 4, and Property 5(b). 
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5. For all possible configurations of covers for the next degree: 

(1) Include the next level to T' to a new poset T", 

(2) Call EnumeratePosets2(r',p+ 1). 

Remark 4.6. One advantage of Algorithm 4.5 is that it is independent of the struc- 
ture of H4. Indeed, we have tested it as well with the root posets of types H3, B4, 
and i*4 (and, of course, adjusted combinatorial information). It turns out that (up 
to symmetries) these posets are uniquely determined by these properties. 

4.2.4. Results. A slight modification of Algorithm 4.5 produces all posets with 
Properties 1 through 4 and satisfying as well Property 5(a). There are 4 such 
posets (see Figure 5), none of them satisfies Property 5(b). We do indeed think 
that the Panyushev action does simply not behave well in type H4. 



4.3. Assuming as well Property 6. In Section 5, we discuss the computations 
that have led to a conjectured Hilbert series of the diagonal coinvariants which 
should possibly provide the numbers of order ideals of given sizes (Property 6). A 
slight modification of Algorithm 4.5 is capable of enumerating all posets with Prop- 
erties 1 through 4 and Property 6. Notice that Property 6 is such a strong condition 
that the computation just takes a few seconds this time. It turns out that there is 
no poset agreeing with the polynomial T-LHi{q,t) conjectured in Conjecture 5.3. 



4.3.1. Assuming Conjecture 5.3 to be false. Even though we strongly believe the 
Hilbert series in Conjecture 5.3 to be true, one might ask how the situation looks like 
if this conjecture turns out to be false. In particular, if it could then be possible 
to find a poset satisfying Properties I through 4 and Property 6. To this end, 
assume now that there exists such a poset, and that the Hilbert series Hff 4 (g,t) 
satisfies assumptions (1) and (2), and that ai. n ^ implies that g*[n] g 2 appears 
in the decomposition of q^HH^iq, q 1 ) ( see Section 5.4 for the definition of the 
two assumptions, of the decomposition of q 60 HH 4 (?, <7~ 1 )j an d of a^.n). By the well 
behavedness with respect to standard parabolic subgroups in Property 1, the terms 
of low degree of Hh 4 (<?, 1) are given by 

I + 4,? + 6q 2 + 7q 3 + 8q A + 8q 5 + 9q 6 + 8q 7 + 8q S + 8q 9 + 9q 10 + ..., 

and the terms of high degree are given by 

... + 2q 49 + q 50 + q 51 + ... + q 60 . 

Then one can compute that there are 180 possible polynomials for {q, t) with 
these properties, which moreover all have different specializations t — 1. Finally, 
there are 120 posets satisfying Properties 1 through 4 and whose order ideals 
yield one of these specializations. Only 7 of the 180 different possible polyno- 
mials T-lHiiq^t) occur. The following list gives the 7 polynomials and the number 
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of corresponding posets. 



2 


X 


[61],,* 


+ g 1 i 1 [49] (?)t - 


hg 3 t 3 [4l] 9 , t 


+ q 1 t 1 [37] q , t 


+ g 4 i 4 [31] 9 , 4 + 








+ g 1 i 1 [21] ?>t - 


hg 2 t 2 [13] g , t 


+ q 6 t 6 


+ q 10 t 1Q 


10 


X 


[61],,t 


+ 9 1 t 1 [49] g ,t- 




+ q 3 t 3 [37] q , t 


+ g 1 * 1 [31] 9 , t + 






q 4 t 4 [25] q:t 


+ g 2 i 2 [21]< M - 


h<? 2 i 2 [13] g , 4 


+ q 6 t 6 


+ q 1Q t 1Q 


12 


X 


[61],,* 


1 1 r . s^-i 

+ q 1 i l [A9} q>t - 


1 1 r j - i 


+ q 4 t i [37} q , t 


+ g 1 t 1 [31] ? , t + 








+ q 2 t 2 [21] qit - 


hg 2 i 2 [13] g ,t 


+ q 6 t 6 


+ q 10 t 10 


1G 


X 


[61],,t 


+ q 1 t 1 [A9] qit - 




+ <Z 4 i 4 [37] g , t 


+ g 1 t 1 [31] ? , t + 






q 2 £ 2 [25] 9:f 


+ q 2 t 2 [21} q , t - 


Hg 3 t 3 [13] g , t 


+ q 6 t 6 


+ q 1Q t 1Q 


20 


X 


[61] fl ,t 


+ q 1 t 1 [A% tt ~ 


Hg 3 i 3 [4l] g , t 


+ q 1 t 1 [37] q , t 


+ « 1 t 1 [31],, t + 






<Z 4 i 4 [25] 9 , t 


+ q 2 t 2 [21] qit - 


hg 2 * 2 [l3] 9)t 


+ q 6 t 6 


+ g 10 i 10 


20 


X 


[61],,t 


+ q 1 t 1 [49] q>t - 




+ q 3 t 3 [37] q j. 


+ < 7 1 t 1 [31] ? , t + 






g 2 t 2 [25] 9 , t 


+ q 2 t 2 [21] q!t - 


H<Z 4 * 4 [13],, t 


+ q 6 t 6 


+ g 10 i 10 


40 


X 


[61kt 


+ q 1 t 1 [49] qtt - 


Hg 3 t 3 [4i] g> t 


+ g 1 * 1 [37] 9 , t 


+ g 1 t 1 [31] g , t + 






g 2 i 2 [25] ? , t 


+ q 2 t 2 [2l] qtt - 


Hg 4 * 4 [i3] 8l i 


+ g 6 * 6 


+ g 10 i 10 



5. Computation of the Hilbert series of the diagonal coinvariants 

To compute the Hilbert series 'M-H i {q,t) = , H{M.{H±);q,i) in Property 6, we 
must compute a minimal basis B of the ideal in C[V © V] generated by the elements 

9{m) := det(ui) w(m) 

for all monomial m = x" 1 ■ ■ ■x% A y{ 1 ■ ■ ■ yf 4 G C[x,y]. Here, W = VF(ff 4 ) is the 
Coxeter group of type H4 acting diagonally on C[x, y]. Each element / G B is 
homogeneous in x and in y of bidegree (d x , d y ), and the Hilbert series is then given 

by 

^(g,t) = $> d ^ dy - 

In this process, we face the following two problems. First, there is the problem 
of computing 0(m) for a given monomial m G C[x, y]. This computation is closely 
related to computing the Reynolds operator in classical computational invariant 
theory (over polynomial rings) since W has a subgroup of index 2 in which all 
elements have determinant 1. However, computing invariants for a group of size 
7200 acting on an 8-dimensional space is slightly out of the range of the classical 
algorithms. Second, there is the problem that computations in fields of coefficients 
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that can be used to provide faithful 4-dimensional turn out to be too slow for the 
needed Grobner basis computations. 

To work around these problems, we are forced to resort to some heuristics. We 
present three steps which accelerate the computation considerably, two of them are 
assumptions under which we cannot prove anymore that the computed series is 
indeed correct. 

5.1. Finding a nice representation. The largest abelian subgroup U in W has 
50 elements. We choose a representation of W in which these 50 elements are 
diagonal matrices. Luckily, for this choice of a basis, W has a subgroup U' of 400 
elements which are all monomial matrices. Thus we compute 0(m) in the following 
way. 

Algorithm 5.1. Theta(m) 
Apply the map 9. 

Input: A monomial m 6 C[xi, . . . , X4, y\, . . . ,2/4]. 

Output: 9(m). 

1. Compute z := J2 w eu det(iy) w(m). This is very easy, since all elements of 
U are diagonal matrices. 

2. If z ^ 0, then: 

(1) Compute z' :— J2 w ueU' /u det(ui) w(m). Computing each of the sum- 
mands is just evaluation, since all elements of U' are monomial matri- 
ces. 

(2) If z' ^ 0, then return g := J2 wU 'ew/U' det(to) w(z'). 

3. Return 0. 

Remark 5.2. The first sum has 50 summands, the second one 8, and the last sum 
has 36 summands. Thus instead of computing 14400 times w(m), we only need 
36 of these expensive computations, and we only perform them if z and z' are not 
zero. In practice, before returning g we also evaluate g at (y\, . . . , 7/4, xx, ■ ■ ■ , X4) to 
avoid a second computation of this invariant. 

5.2. Choosing a fast field of coefficients. In principle, faithful 4-dimensional 
representations of W require a field containing the golden ratio, for instance Q(\/5). 
To apply Algorithm 5.1, we need fifth roots of unity because we need to diagonalize 
simultaneously the subgroup U. This is not a problem since implementations of 
cyclotomic fields are highly optimized in most computer algebra systems. However, 
although they are quite fast, these fields (in characteristic zero) are too slow for 
the last step, the computation of the Grobner basis. We therefore perform all 
computations over a finite field ¥ q of order q with 5 | (q — 1) (so that the 5-th roots 
of unity are included). This can result in a slightly different Hilbert series, but with 
large q this is very unlikely. 

5.3. Reducing the number of variables. Finally, since we are only interested in 
the degrees of the components of the elements of B, we may evaluate y\, . . . ,1/4 at 
Z\t, . . . , Zit for some zi, . . . , £4 G ¥ q before computing B. This could lead to some 
wrong coefficients in "Hff 4 (g,t). But without this simplification, the computation 
of B is not feasible on modern computers. 
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5.4. Further assumptions on rl LM ($); q, t ) . Even though we already have made 
concessions, the above techniques are not quite sufficient to obtain all the required 
coefficients. We finally fill the gaps in the computed series by using the following 
further assumptions. 

(1) The series H(M($);q,t) is of the form 

a itn qH l [n] qt t 

ieN,neN 

for some a^„ G N, i, n e N. 

(2) The evaluation t = q^ 1 in the Hilbert series yields 2 

rfi*U (« - [32] g [42] g [50] g [60] g 
q n HM ,q ) [2]g[12]q[20]q[mq , 

where we write [n] q = [n\ q ,i = q" -1 + q' l ~ 2 + . . . + 2 + 1. 

Under these assumptions, the above heuristics, and with extensive computations 
we obtain the following conjecture. 

Conjecture 5.3. The Hilbert series of the diagonal coinvariants of the Coxeter 
group of type H4 is given by 

H Ht (q,t) = [61} q , t + qt[49] q , f + «t[41] g , t + q 2 t 2 [37] qA + ?t[31] git + 
qH 3 [25] q<t qH 2 [2l] qtt + qH^li]^ + q e t 6 + q W t W . 

Notice that the polynomial q 60 'Hii i (q, q -1 ) in assumption (2) has the decompo- 
sition 

q 60 n H Aq,q^) = [61],* +<7 12 [49] g2 + g 20 [4l] 9 2 + g 24 [37] g 2 +<7 30 [31] 9 2 + 

< Z 36 [25] g2 + (? 40 [21] 92 + (7 48 [13] g2 +2g 60 , 

and that this decomposition in such summands is unique in the sense that consec- 
utive summands g Q [6] g 2 and <7 c [d] g 2 satisfy a < c and b > d. One can now observe 
that the sequence 

61,49,41,37,31,25,21,13,1,1 

of numbers involved in this decomposition coincides with the sequence of numbers 
involved in the conjectured Hilbert series. 
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